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SECTION  1 
INTRODUCTION 


There  are  currently  many  different  computer  codes  which  can  be  used  for 
computations  of  intense  impulsive  loading  due  to  high  velocity  impact  and/or 
explosive  detonation.  Although  the  current  status  of  these  codes  is  that  they  can 
now  be  used  to  perform  meaningful  computations,  it  is  generally  agreed  that  there  is 
a  need  for  improved  strength  and  fracture  models.  There  is  also  a  need  to  develop 
efficient  procedures  to  obtain  constants  for  these  models. 

The  work  described  in  this  report  is  focused  on  two  computational  strength 
models:  The  Johnson-Cook  model"  and  the  Zerilli-Armstrong  model.2  The  report 
first  describes  the  models  and  discusses  relative  comparisons.  Then  model 
predictions  are  compared  with  experimental  data  from  tension/torsion  tests  and 
cylinder  impact  tests.  This  is  followed  by  a  discussion  of  how  cylinder  impact  tests 
may  be  used  to  determine  constants  for  these  models.  The  report  ends  with  a 
summary  and  conclusions. 
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SECTION  2 

THE  JOHNSON-COOK  AND  ZERILLI-ARMSTRONG  MODELS 


Most  computational  strength  models  express  the  equivalent  (von  Mises)  tensile 
flow  stress  as  a  function  of  the  equivalent  plastic  strain  rate,  temperature  and/or 
pressure.  For  the  Johnson-Cook  model1  the  equivalent  tensile  flow  stress  is 
expressed  as 

o  =  [A  +  Bcn][l  +  Cine*]  11  —  T*m]  (1) 

where  e  is  the  equivalent  plastic  strain,  i*  =  c/eq  is  the  dimensionless  plastic  strain 
rate  for  C0  =  1.0s'1,  and  T*  =  (T-TroomV(Tmelt-^oom)  is  the  homologous 
temperature.  The  relationship  is  valid  for  0  ^  T*  <  1.0.  The  five  material  constants 
are  A,  B,  n,  C,  and  m. 

The  expression  in  the  first  set  of  brackets  gives  the  stress  as  a  function  of  strain 
for  t*  =  1.0  and  T*  =  0.  The  expressions  in  the  second  and  third  sets  of  brackets 
represent  the  effects  of  strain  rate  and  temperature,  respectively.  At  the  melting 
temperature  (T*  =  1.0),  the  stress  goes  to  zero  for  all  strains  and  strain  rates.  The 
basic  form  of  the  model  is  readily  adaptable  to  most  computer  codes,  since  it  uses 
variables  (e,  i *,  T*)  which  are  available  in  the  codes. 

This  model  has  some  desirable  features  inasmuch  as  it  is  simple  to  implement, 
does  not  require  excessive  computing  time  or  memory,  can  be  used  for  a  variety  of 
metals,  constants  can  be  straightforwardly  obtained  from  a  limited  number  of 
laboratory  tests,  and  the  effects  of  the  important  variables  are  readily  identifiable 
and  separable. 

The  primary  disadvantage  of  this  model  is  that  it  is  empirical  and,  therefore, 
has  no  sound  physical  basis.  This  means  that  exceptional  care  must  be  exercised 
when  using  it  for  extrapolated  values  of  e,  £*,  and  T*. 

Some  of  the  motivation  and  background  for  this  model  can  be  obtained  from 
References  3  to  6,  where  an  attempt  was  made  to  understand  and  simulate  tests  of 
large  torsional  strains  over  a  range  of  strain  rates. 

Although  various  test  techniques  can  be  used  to  obtain  constants  for  this 
model,  the  following  has  worked  well.1  First,  the  yield  and  strain  hardening  con¬ 
stants  (A,  B,  n)  are  obtained  from  isothermal  tension  and  torsion  tests  at  relatively 
low  strain  rates  (i*  <  1.0).  Next,  the  strain  rate  constant,  C,  is  determined  from 
torsion  tests  at  various  strain  rates,  and  tension  tests  (quasi-static  and  Hopkinson 
bar)  at  two  strain  rates.  Finally,  the  thermal  softening  constant,  m,  is  determined 
from  Hopkinson  bar  tests  at  various  temperatures. 
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Sometimes  it  is  possible  to  obtain  the  yield  and  strain  hardening  constants 
(A,  B,  n)  from  cylinder  impact  tests.  This  will  be  more  fully  discussed  in  Section  4. 

The  Zerilli- Armstrong  model  is  based  on  dislocation  mechanics,  and  is,  there¬ 
fore,  more  physically  based  than  the  Johnson-Cook  Model.  The  Zerilli- Armstrong 
model  has  two  forms:  one  for  face  centered  cubic  (fee)  metals;  and  another  for  body 
centered  cubic  (bcc)  metals.2  The  expression  for  fee  metals  is 

o  =  C0  +  C2e1/2exp(-C3T  +  C4Tlne)  (2) 

where  e  is  the  equivalent  plastic  strain,  £  is  the  equivalent  strain  rate,  and  T  is  the 
absolute  temperature.  The  four  constants  are  Cq,  C2,  C3,  and  C..  Here  the  initial 
yield  stress,  C0,  is  independent  of  strain  rate  ana  temperature.  Also,  the  stress  does 
not  necessarily  go  to  zero  at  the  melting  temperature.  Reference  2  provides  a 
discussion  of  how  CQ  is  affected  by  solute  and  grain  size. 

The  expression  for  bcc  metals  is 

o  =  C0  +  Cj  exp  (-  C3T  +  C4T  In  c}  +  C&en  (3) 

where  the  variables  (e,  £,  T)  are  as  defined  for  Equation  (2)  and  six  constants  are  CQ, 
C,,  Co,  C4,  C5,  and  n.  Here  the  initial  yield  stress  is  a  function  of  CQ,  C. ,  C3,  and  C4. 
Again,  the  stress  does  not  necessarily  go  to  zero  at  the  melting  temperature. 

It  can  be  seen  in  Equation  (2)  that  the  strain,  strain  rate,  and  temperature 
effects  are  all  coupled  together  for  the  fee  model.  In  Equation  (3),  however,  the  effect 
of  strain  hardening  is  separated  from  the  coupled  strain  rate  and  temperature  for  the 
bcc  model. 

The  fact  that  the  Zerilli- Armstrong  model  is  based  on  dislocation  mechanics 
makes  it  preferable  to  the  Johnson-Cook  model.  On  the  other  hand,  the  more  complex 
form  of  tiie  Zerilli-Annstrong  model  appears  to  make  it  more  difficult  to  obtain  the 
appropriate  constants. 

Figure  1  shows  a  comparison  of  isothermal  and  adiabatic  stress-strain  relation¬ 
ships  for  the  two  models  of  interest.  This  is  done  for  OFHC  Copper  (fee)  and  Armco 
Iron  (bcc).  Constants  for  both  models  were  obtained  from  essentially  the  same  data 
base.1*2 

For  the  OFHC  copper,  the  Johnson-Cook  model  predicts  higher  adiabatic 

stresses  at  lower  strains,  and  lower  stresses  at  higher  strains.  For  both  the  OFHC 
copper  and  the  Armco  iron,  the  Johnson-Cook  model  predicts  less  of  a  strain  rate 
effect  than  does  the  Zerilli- Armstrong  model.  Generally,  however,  the  adiabatic 
responses  are  similar  for  the  two  models. 

Figure  2  shows  the  effect  of  strain  rate  for  the  two  models.  For  the  Armco  Iron 
the  Zerilli- Armstrong  model  shows  a  much  stronger  reliance  on  strain  rate  at  the 
higher  strain  rates  shown.  This  is  consistent  with  torsion  data  for  Armco  IF  Iron5 
and  with  trends  predicted  by  other  researchers. 
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FIGURE  1.  ISOTHERMAL  AND  ADIABATIC  STRESS-STRAIN  RELATIONSHIPS  FOR 
OFHC  COPPER  AND  ARMCO  IRON  USING  THE  JOHNSON-COOK 
AND  ZERILLI-ARMSTRONG  MODELS 
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FIGURE  2.  ADIABATIC  STRESS-STRAIN  RATE  RELATIONSHIPS  FOR  OFHC 
COPPER  AND  ARMCO  IRON  USING  THE  JOHNSON-COOK 
AND  ZERILLI-ARMSTRONG  MODELS 
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SECTION  3 

COMPARISON  OF  MODEL  PREDICTIONS  AND  TEST  DATA 


An  assessment  of  the  models  can  be  made  by  comparing  the  model  predictions  to 
experimental  data.  Figure  3  shows  four  comparisons  for  OFHC  copper.  The  upper 
left  portion  of  Figure  3  is  for  dynamic  Hopkinson  bar  test  data  at  room  temperature, 
and  the  upper  right  portion  is  for  similar  data  at  an  elevated  temperature.  It  is 
assumed  that  the  responses  are  adiabatic.  There  is  relatively  good  agreement 
between  the  model  predictions  and  the  test  data. 

The  lower  portion  of  Figure  3  shows  comparisons  with  quasi-static  tension  data 
and  quasi-static  torsion  data.  Here  the  responses  are  assumed  to  be  isothermal.  Both 
models  tend  to  underpredict  the  strength  in  tension  and  overpredict  the  strength  in 
torsion.  This  is  due,  in  part  to  the  fact  that  real  materials  do  not  always  obey  the  von 
Mises  flow  rule,5  which  is  an  inherent  part  of  most  hydrocode  computational 
algorithms.  The  von  Mises  flow  rule  states  that  the  equivalent  tensile  stress  and 
strain  are  a  =  V3x  and  e  =  y/V3,  where  x  and  y  are  the  shear  stress  and  strain. 

Figure  4  shows  similar  comparisons  for  Armco  iron.  Here  again  there  is  good 
general  agreement  with  the  Hopkinson  bar  data.  Both  models  tend  to  underpredict 
the  strength  for  the  quasi-static  tension  and  torsion  data.  The  reason  this  occurs  for 
the  Johnson-Cook  model  is  that  the  primary  application  is  for  higher  strain  rates, 
and  the  strain  rate  constant  was  therefore  selected  to  give  better  correlation  with  the 
higher  strain  rate  torsion  data.  This  results  in  an  underprediction  of  strength  for 
lower  strain  rates  (i*  <  1.0). 

Another  interesting  evaluation  is  to  compare  the  computed  shapes  of  cylinder 
impacts  onto  a  rigid  surface,  with  actual  test  data,  as  shown  in  Figure  5.  To  quantify 
the  degree  of  agreement  between  computed  shapes  and  test  data,  an  average  error 
has  been  defined  as 


|ALj  jAD|  |AW| 

IT  +  D  +  w 


(4) 


where  L,  D,  and  W  are  the  deformed  length,  diameter,  and  bulge  (diameter  at  0.2Lq 
from  the  deformed  end)  from  the  test  results,  and  AL,  AD,  and  AW  are  the  differences 
between  the  computed  and  test  results.  It  can  be  seen  that  both  models  give  good 
general  agreement,  but  that  the  Zerilli- Armstrong  model  gives  better  agreement.2 

The  maximum  computed  strains  in  Figure  5  fall  within  the  range 
1.57  <  e  x  s  2.04.  The  strain  distributions  in  Figure  6,  however,  show  that  the 
overwhelming  majority  of  the  elements  experience  equivalent  strains  less  than  0.6. 
Therefore,  the  comparisons  in  Figure  5  tend  to  reflect  the  accuracy  of  a  model  for 
relatively  low  strains  (e  <  0.6),  but  do  not  necessarily  provide  a  good  indication  for 
larger  strains  (e  >  0.6). 
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FIGURE  3.  COMPARISON  OF  JOHNSON-COOK  AND  ZERILLI-ARMSTRONG  MODEL 
RESULTS  WITH  OFHC  COPPER  TEST  DATA 
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FIGURE  4.  COMPARISON  OF  JOHNSON-COOK  AND  ZERILLI-ARMSTRON G  MODEL 
RESULTS  WITH  ARMCO  IRON  TEST  DATA 
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FIGURE  5.  COMPARISON  OF  CYLINDER  IMPACT  TEST  RESULTS  AND  COMPUTED 
SHAPES  FOR  OFHC  COPPER  AND  ARMCO  IRON,  USING  THE 
JOHNSON-COOK  AND  ZERILLI- ARM STRON G  MODELS 
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FIGURE  6.  DISTRIBUTION  OF  COMPUTED  STRAINS  OF  OFHC  COPPER  AND 
ARMCO  IRON  CYLINDER  IMPACT  COMPUTATIONS,  USING 
THE  JOHNSON-COOK  AND  ZERI LLI-ARMSTRON  G  MODELS 
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SECTION  4 

DETERMINATION  OF  MODEL  CONSTANTS  FROM  CYLINDER 

IMPACT  TEST  DATA 


The  preceding  sections  have  attempted  to  define,  understand,  and  evaluate  the 
Johnson-Cook  and  Zerilli-Armstrong  strength  models.  The  constants  for  these 
models  were  obtained  primarily  from  torsion  tests  at  various  strains  rates, 
Hopkinson  bar  tests  at  various  temperatures,  and  quasi-static  tension  tests. 
Performing  these  tests  can  often  be  time-consuming  and/or  expensive. 

The  cylinder  impact  test,  on  the  other  hand,  is  simple,  inexpensive,  and 
exhibits  large  strains,  high  strain  rates,  and  elevated  temperatures.  Unfortunately, 
the  strains,  strain  rates,  temperatures,  and  stresses  vary  throughout  the  test 
specimen  and  throughout  the  duration  of  the  test. 

During  the  past  years,  there  have  been  several  attempts  at  defining  dynamic 
flow  stresses  from  cylinder  impact  test  results.7"10  Figure  7  shows  how  cylinder 
impact  test  data  can  be  used  to  obtain  constants  for  the  Johnson-Cook  Model.  The 
test  data  are  identical  to  those  shown  in  Figure  5.  The  adiabatic  stress-strain 
relationships  are  shown  in  Figure  8  and  the  strength  constants  are  given  in  Table  1. 
The  adiabatic  stresses  are  shown  only  to  the  maximum  strains  attained  in  the 
computed  results.  Also  included  in  Table  1  and  Figure  8  are  data  from  Reference  1 
and  Figure  1. 

For  Case  A-l  the  length,  L,  of  the  deformed  cylinder  (OFHC  copper)  is  matched 
with  the  computational  result.  Because  only  one  deformed  dimension  is  matched, 
only  one  independent  strength  constant  (representing  a  constant  flow  stress)  can  be 
obtained.  It  can  be  seen  that  there  are  significant  discrepancies  between  the  test  and 
computational  results  at  the  deformed  end  of  the  cylinder. 

Case  B-l  allows  for  linear  strain  hardening.  By  matching  both  the  deformed 
length,  L,  and  the  maximum  diameter,  D,  it  is  possible  to  obtain  the  two  constants 
for  the  linear  hardening.  Here  the  computed  result  is  in  good  general  agreement 
with  the  test  result.  If  this  model  would  be  applied  to  larger  strains  than  those 
experienced  in  the  computed  results  of  Case  B-l,  however,  it  would  probably  lead  to 
excessively  high  stresses.  It  is  well  known  that  thermal  softening  tends  to  decrease 
the  rate  of  strain  hardening  at  large  strains,  as  shown  in  Figure  1. 

Cases  C-l  and  D-l  are  based  on  the  Johnson-Cook  model  of  Equation  (1). 
Although  there  are  five  constants  in  this  model  (A,  B,  n,  C,  m),  only  the  yield  and 
strain  hardening  constants  (A,  B,  n)  are  determined  from  the  test  data.  The  strain 
rate  constant,  C,  and  the  thermal  softening  constant,  m,  must  be  approximated  or 
obtained  from  other  sources.  Because  there  are  now  three  constants,  it  is  possible  to 
match  three  deformed  dimensions:  the  length,  L,  the  maximum  diameter,  D,  and  an 
intermediate  diameter,  W. 
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OFkC  COPPER 

i - 1 

i  i 

CASE  A-1 


I - T  i - 1  l - 1 

i  i  II  li 

CASE  0-1  CASE  C-1  CASE  D-1 


NOTES: 

•  OFHC  COPPER  CYLINDER  (L„  =  25.4mm,  D0  »  7.9nm,  LTE8T  =  10.2mm,  DTEST  =  13.5mm,  WTEST  =  10.1mm) 
•ARMCO  RON  CYLINDER  (L q  —  25.4mm,  Do  =  2.6mm,  L^g^s  10.8mm,  Dtests  13.7mm,  W TE3T  ~  8.6mm) 

•  'WPACT  VELOC1TES  =  190  m/s  (OFHC  COPPER)  AND  221  m/s  (ARMCO  RON) 

•  TEST  RESULTS  HDIGATED  BY  DOTS  t  •  » 

•  EQUIVALENT  PLASTIC  STRAW  CONTOURS  SHOWN  AT  .1,  E,  .4.  .5  6  1X1,  1.2,  1.5 

•  ALL  CASES  COMPUTED  WfTH  JOHNSON  -  COOK  MODEL 


FIGURE  7.  EXAMPLES  OF  CYLINDER  IMPACT  TEST  RESULTS  USED  TO  DETERMINE 
OFHC  COPPER  AND  ARMCO  IRON  CONSTANTS  FOR  THE  JOHNSON- 
COOK  MODEL 


equivalent  flow  STRESS,  a  (MPa) 


NSWCTR  88-250 


O  0.5  1.0  1.5  2.0  2j5  3.0  4.5  5X)  5.5 


EQUIVALENT  PLASTIC  STRAW,  E 


FIGURE  8.  COMPARISON  OF  ADIABATIC  STRESS-STRAIN  RELATIONSHIPS  FOR 
VARIOUS  OFHC  COPPER  AND  ARMCO  IRON  CONSTANTS , 

USING  THE  JOHNSON-COOK  MODEL 
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TABLE  1.  OPHC  COPPER  AND  ARMCO  IRON  CONSTANTS  FOR  THE  JOHNSON-COOK 
MODEL,  AS  DETERMINED  FROM  CYLINDER  IMPACT  TEST  RESULTS 


CONSTANTS  FOR  JOHNSON  -  COOK  MODEL 
-  [A+B£n  ]  [l+C!nS*  J  [l-T*  m] 

OFHC  COPPER 

CASE  A-1 

CASE  B-1 

CASE  C-1 

CASE  D-1 

REF.  1 

A 

(Mpa) 

250 

mm 

98 

90 

B 

(Mpa) 

0. 

- 

368 

292 

n 

-  • 

.70 

.31 

C 

0. 

mgm 

.025 

.025 

m 

- 

■■ 

1.00 

1.09 

1.09 

ARMCO  IRON 

CASE  A-2 

CASE  B-2 

CASE  C-2 

CASE  D-2 

REF.1 

A 

(Mpa) 

593 

NOT  POSSIBLE  TO 

175 

B 

(Mpa) 

0. 

OBTAIN  CONSTANTS 

380 

n 

- 

1.0 

FOR  CASES  C-2 

.32 

C 

0. 

0. 

AND  D-2 

.060 

m 

•  - 

- 

.55 

S*-!/e0  for  !o 

=  1.0  S'1 

Ij  T  «*(T-T room ) / (Tmelt  -T room  ) 
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It  should  be  noted  that  D  and  W  are  not  totally  independent  of  one  another.  For 
a  given  length,  L,  W  will  tend  to  decrease  as  D  increases.  This  is  simply  due  to  an 
approximate  conservation  of  the  volume  of  the  cylinder. 

Case  C-l  is  the  result  obtained  if  nothing  is  known  about  the  strain  rate  or 
thermal  softening  characteristics  of  the  material.  Here,  the  strain  rate  constant  is 
set  to  C  =  0  and  the  thermal  softening  is  assumed  to  be  linear,  or  m  =  1.0.  Case  D-l 
uses  the  previously  determined  strain  rate  and  thermal  softening  constants  from 
Reference  1  and  Figure  1.  Even  though  there  is  little  apparent  difference  between 
Cases  C-l  and  D-l  for  this  problem,  for  some  other  problems  (with  wider  ranges  of 
strain  rates  and  temperatures)  the  Case  D-l  constants  will  probably  provide  better 
results. 

Looking  at  Figure  8,  it  can  be  seen  that  there  are  significant  discrepancies  in 
the  various  adiabatic  stress-strain  relationships,  especially  at  the  larger  strains. 

Yet,  with  the  exception  of  Case  A-l,  they  all  provide  a  generally  good  correlation  with 
the  test  data.  The  reason  for  this  situation  is  provided  in  the  distribution  of  strain,  as 
shown  previously  in  the  upper-left  portion  of  Figure  6.  Although  some  of  the 
elements  experience  an  equivalent  plastic  strain  as  high  as  2.04,  most  of  the  elements 
experience  strains  less  than  0.6.  Looking  back  to  this  specific  range  of  strains  in 
Figure  8,  there  is  good  general  agreement  (except  for  Case  A-l)  between  the  various 
adiabatic  stress-strain  relationships.  Therefore,  because  the  various  models 
essentially  agree  with  one  another  in  this  relatively  narrow  band  of  strain,  it  would 
be  expected  that  they  would  give  similar  results  for  computed  solutions  whose  strains 
fall  within  this  narrow  band. 

In  the  lower  portion  of  Figure  7,  the  same  approach  was  attempted  for  the 
Armco  iron.  For  Case  A-2  only  the  length,  L,  was  matched,  giving  a  constant  flow 
stress.  For  Case  B-2  the  length,  L,  and  diameter,  D,  were  matched  using  linear  strain 
hardening.  The  resulting  adiabatic  flow  stresses  are  shown  in  Figure  8.  There  is 
clearly  less  strain  hardening  in  the  Armco  iron  than  in  the  OFHC  copper. 

For  Cases  C-2  and  D-2,  however,  it  was  not  possible  to  match  all  three 
dimensions  (L,  D,  W)  by  varying  the  yield  and  strain  hardening  constants  (A,  B,  n). 
This  is  probably  due  to  the  fact  that  there  is  less  strain  hardening  in  the  Armco  iron 
(when  compared  to  the  OFHC  copper),  and  therefore  less  of  a  bulge,  W.  It  appears 
that  it  simply  is  not  possible  to  extract  strain-hardening  constants  for  test  data  which 
exhibit  veiy  little  strain  hardening.  Also,  the  empirical  fonn  of  Equation  (1)  may  not 
be  well  suited  to  accurately  represent  the  Armco  iron. 

For  the  Johnson-Cook  model,  it  appears  that  it  is  possible  to  obtain  the  yield 
and  strain  hardening  constants  (A,  B,  n)  for  some  materials,  but  not  for  all  materials. 
This  is  probably  dependent  on  the  degree  of  strain  hardening  in  the  material.  The 
Strain  rate  and  thermal  softening  constants  (C,  m)  must  be  approximated  or  obtained 
from  other  sources.  It  is  also  not  desirable  to  extrapolate  the  use  of  the  model  to 
larger  strains  than  experienced  in  the  cylinder  impact  test. 

Looking  at  the  Zerilli- Armstrong  models  in  Equations  (2)  and  (3),  it  can  be  seen 
that  four  constants  are  needed  for  fee  metals  (such  as  OFHC  copper)  and  six  con¬ 
stants  are  needed  for  bcc  metals  (such  as  Armco  Iron).  Based  on  previous  experience 
with  the  Johnson-Cook  model,  only  two  or  three  constants  can  be  obtained  from  the 
shape  of  the  deformed  cylinder. 
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Looking  again  at  Equation  (2)  for  fee  metals,  it  would  appear  that  the  first  two 
constants  (C0,  C  J  could  be  obtained  from  the  cylinder  impact  test  data,  if  the  strain 
rate  and  thermal  softening  constants  (C3,  C4)  could  be  approximated  or  obtained  from 
other  sources.  These  results  are  shown  as  Case  E  in  Figure  9.  The  length,  L,  and 
diameter,  D,  have  been  matched  exactly,  and  the  resulting  bulge,  W,  is  very  close  to 
that  of  the  test  data.  The  new  constants,  (C„  and  C2)  were  obtained  in  conjunction 
with  the  previous  values  of  C3  and  C4,  as  taken  from  Reference  2. 

Another  approach  is  to  obtain  the  yield  stress,  C0,  from  another  source,  and  then 
solve  for  C2,  C,,  and  C4  by  matching  the  three  deformed  dimensions  (L,  D,  W).  The 
yield  stress  in  Equation  (2)  is  independent  of  both  the  strain  rate  and  the 
temperature,  and  can,  therefore,  be  easily  determined  from  simple  tension  tests 
and/or  handbook  data.  The  results  of  this  approach  are  shown  as  Case  F  in  Figure  9, 
where  the  new  constants  (C2,  C3,  C4)  were  obtained  for  the  previous  value  of  CQ,  as 
taken  from  Reference  2.  A  summary  of  the  constants  for  Cases  E  and  F  is  given  in 
Table  2. 

It  is  interesting  to  compare  the  adiabatic  stress-strain  relationships  of  Case  D-l 
(Johnson-Cook)  and  Case  F  (Zeri  lli- Armstrong),  because  both  of  these  cases  give  a 
perfect  fit  (A  =  0)  to  the  test  data.  These  two  cases,  along  with  three  other  cases,  are 
shown  in  Figure  10.  As  expected,  the  responses  of  Cases  D-l  and  F  are  very  similar  at 
smaller  strains  (e  <  0.6).  They  are  also  in  good  general  agreement  at  larger  strains. 
Case  E,  where  only  two  constants  were  determined  for  the  Zeri  lli -Armstrong  model, 
is  also  shown  in  Figure  10. 

It  should  be  noted  that  the  actual  values  of  the  constants  in  Cases  E  and  F  may 
or  may  not  be  physically  realistic.  For  these  cases,  the  constants  are  selected  only  to 
match  the  deformed  shape  of  the  cylinder.  The  other  two  curves  are  from  Reference  1, 
which  presents  the  initial  Johnson-Cook  model  and  data,  and  from  Reference  2, 
which  presents  the  initial  Zerilli -Armstrong  model  and  data. 

An  attempt  was  made  to  determine  three  constants  (Cft,  C.,  n)  from  the  Zerilli- 
Armstrong  model  of  Equation  (3)  for  Armco  iron,  using  the  reported  values2  of  the 
other  three  constants  (C2,  C3,  C4).  As  was  the  case  with  the  Johnson-Cook  model 
(Cases  C-2  and  D-2),  however,  it  was  not  possible  to  match  all  three  dimensions 
(L,  D,  W)  of  the  test  results. 
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l - 1 

i  I 

CASE  E 


I - I 

!  I 

CASE  F 


A  =  .001 
* 

6maxs1-59 


NOTES: 

•  OFHC  COPPER  CYLINDER  <L0  =  25.4mm,  D0  =  7.6mm, 

Ltest  =  16-2mm-  dtest  =  13-5mm  .  wiEST  = 

•  IMPACT  VELOCITY  =  190  M/S 

•  TEST  RESULTS  INDICATED  BY  DOTS  •  •  • 

•  EQUIVALENT  PLASTIC  STRAIN  CONTOURS  SHOWN  AT 
.1,  .3,  .4,  .5,  .6,  1.0,  1.2,  1.5 

•  BOTH  CASES  COMPUTED  WITH  ZERILLhARMSTRONG  MODEL 


FIGURE  9.  EXAMPLES  OF  CYLINDER  IMPACT  TEST  RESULTS  USED  TO  DETERMINE 
OFHC  COPPER  CONSTANTS  FOR  THE  ZERI LLI- ARM STRON G  MODEL 


19 


NSWCTR  88-250 


TABLE  2.  OFHC  COPPER  CONSTANTS  FOR  THE  ZERILLI-ARMSTRONG 
MODEL,  AS  DETERMINED  FROM  C YLIN  DER  IMPACT 
TEST  RESULTS 


CONSTANTS  FOR  ZERILLI-ARMSTRONG 
FArp  rPMTcncn  r.i  ibio.  unnp 

OFHC  COPPER 

G»Co+C2  8  0,5 

EXP(-C3T+C4T  Am  8) 

CASEE 

CASE  F 

REF  2 

Co 

(Mpa)- 

38 

65 

65 

C2 

(Mpa) 

903 

546 

890 

C3 

(»K-1  ) 

.0028 

.0025 

.0028 

C4 

(oK-1  ) 

.000115 

.000199 

.000115 
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EQUIVALENT  FLOW  STRESS,  <i  (Mpa) 
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FIGURE  10.  COMPARISON  OF  ADIABATIC  STRESS-STRAIN  RELATIONSHIPS 
FOR  VARIOUS  OFHC  COPPER  CONSTANTS,  USING  BOTH  THE 
JOHNSON-COOK  AND  ZERILLI-ARMSTRONG  MODELS 
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SECTION  5 

SUMMARY  AND  CONCLUSIONS 


This  report  has  attempted  to  describe  and  evaluate  two  computational  strength 
models,  and  to  determine  how  constants  might  be  obtained  from  cylinder  impact  test 
results.  Some  conclusions  are  as  follows: 

•  The  Zerilli- Armstrong  model  is  more  physically  based  and  provides  better 
agreement  with  the  cylinder  impact  test  data. 

•  The  Johnson-Cook  model  has  a  more  simple  form  which  allows  the 
constants  to  be  obtained  in  a  more  straightforward  manner. 

•  Both  models  show  generally  good  agreement  with  Hopkinson  bar  test  data 
at  relatively  low  strains. 

•  Both  models  show  discrepancies  with  quasi-static  tension  and  torsion  data 
at  large  strains.  This  is  probably  due  to  the  inadequacy  of  the  von  Mises 
flow  rule  at  large  strains. 

•  Cylinder  impact  test  results  can  be  used  to  determine  two  or  three 
constants  for  either  the  Johnson-Cook  model  or  the  Zerilli- Armstrong 
model.  The  other  constants  must  be  estimated  or  obtained  from  other 
sources. 

•  When  the  length  and  diameter  of  the  cylinder  impact  test  data  are 
matched  by  determining  two  constants,  the  resulting  bulge  is  generally  in 
close  agreement  because  the  cylinder  tends  to  conserve  volume. 

•  The  accuracy  of  either  model  has  not  been  evaluated  for  the  combination  of 
large  strains  (e  >  0.6)  and  high-strain  rates  (c  >  103).  In  unpublished 
work  not  performed  in  this  contract,  however,  the  Johnson-Cook  model  has 
been  used  to  give  good  overall  predictions  of  Explosive  Formed  Penetrators 
for  both  OFHC  copper  and  Armco  iron.  Here,  most  of  the  material  is 
strained  to  a  much  greater  extent  than  it  is  in  the  cylinder-impact  tests. 
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